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Abstract 

Time-correlated single photon counting has recently been combined with mode- 
locked picosecond pulsed excitation to measure the fluorescent lifetimes and energy 
emissions of single molecules in a flow stream. Maximum likelihood (ML) and least 
squares methods agree and are optimal when the number of detected photons is large 
however, in single molecule fluorescence experiments the number of detected photons 
can be less than 20, 67% of those can be noise and the detection time is restricted to 1 0 
nanoseconds. Under the assumption that the photon signal and background noise are two 
independent inhomogeneous Poisson processes, we derive the exact joint arrival time 
probability density of the photons collected in a single counting experiment performed in 
the presence of background noise. The model obviates the need to bin experimental data 
for analysis, and makes it possible to analyze formally the effect of background noise on 
the photon detection experiment using both ML or Bayesian methods. For both methods 
we derive the joint and marginal probability densities of the fluorescent lifetime and 
fluorescent emission. The ML and Bayesian methods are compared in an analysis of 
simulated single molecule fluorescence experiments of Rhodamine 1 1 0 using different 
combinations of expected background noise and expected fluorescence emission. While 
both the ML or Bayesian procedures perform well for analyzing fluorescence emissions, 
the Bayesian methods provide more realistic measures of uncertainty in the fluorescent 
lifetimes. The Bayesian methods would be especially useful for measuring uncertainty in 
fluorescent lifetime estimates in current single molecule flow stream experiments where 
the expected fluorescence emission is low. Both the ML and Bayesian algorithms can be 
automated for applications in molecular biology. 
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Introduction 

In the last four years it has become possible to record temporal data from 
individual fluorescent dye molecules in flow streams and to estimate single molecule 
fluorescent lifetime and energy emissions (21). Perfection of this technology should 
significantly enhance the analytic sensitivity of fluoroimmunoassays, capillary zone 
electrophoresis, flow cytometry, DNA fingerprinting, fragment sizing sequencing. Proper 
application of the single molecule fluorescence techniques in these detection systems 
requires accurate assessment of the uncertainty in the measurement of the molecule's 
lifetime and energy emission. As a consequence, there is much interest in characterizing 
the statistical properties of photon data collected in single molecule excitation 
experiments (19). 

Hall and Selinger (5) studied the general problem of decay lifetime estimation and 
derived equations for the maximum likelihood (ML), method of moments and least 
squares (LS) estimations of the parameters in an inhomogeneous Poisson process model 
of photon detection in the absence of background noise. They also reported large sample 
formulae for the variances of these estimators, and showed that the ML estimation 
procedure was the most statistically efficient and that it was also mathematically 
tractable. Peck et al. studied single molecule fluorescence detection with the signal and 
background noise processes modelled as two independent homogeneous Poisson 
processes and described an autocorrelation procedure for burst detection. Tellinghuisen 
and Wilkerson studied the performance of ML and LS methods for estimating 
fluorescent lifetimes using an arrival time model in the case where N, the numbers of 
detected photons, is small and T, the observation time, is finite. They found that both 
estimates were biased and that the variance of the ML estimate was formally divergent. 
They suggest that the (N/ N- 1) bias in the reciprocal of the lifetime parameter could be 
easily corrected and that statistical models parameterized in terms of this parameter were 
more analytically tractable. Tellinghuisen studied the properties of least squares (LS) 
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techniques for estimating fluorescent lifetimes in the presence of background noise using 
Monte Carlo methods and partial derivative matrices to perform error analyses. He 
reported that background noise reduced the precision in the lifetime estimate by an order 
of magnitude and that this measurement precision could be recovered by incorporating 
into the analysis a model of the background noise based on the properties of the detection 
system studied in the absence of a fluorescent experiment. 

In the study of photon bursts from single Rhodamine 1 10 dye molecules dissolved 
in methanol, Tellinghuisen et al. compared the use of an approximate arrival time plus 
background noise model analyzed with LS methods to an arrival time noise free model 
analyzed by ML procedures. In their LS analysis the error in the reciprocal of the lifetime 
parameter exceeded the N* expected for data collected over an infinite time interval in 
the absence of background noise. Because the model analyzed with the LS methods 
included background noise and because problems in the parameter estimation could be 
readily diagnosed with contour plots of the minimum y} function, these authors 
concluded that the LS approach offered the preferred means of analyzing the effect of 
background noise in single molecule fluorescence experiments. 

The principal analytic issues to be addressed are: (1) proper specification of a 
probability model for photon detection from single molecules in the presence of 
background noise; and (2) formulation of a statistical estimation procedure based on that 
model which allows accurate determination of single molecule fluorescent emissions and 
lifetime when emissions are low background noise is significant and the detection 
interval is finite. There is significant evidence to suggest that detection of photon from 
single fluorescent molecule may be modelled a Poisson process. Therefore, under the 
assumption that the signal and background noise are two independent inhomogeneous 
Poisson processes, we derive the exact joint probability density of photon arrivals in a 
finite time interval from a single fluorescent molecule derive both ML and Bayesian 
procedures for estimating the fluorescence lifetime energy emission and their associated 
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uncertainties in terms of well defined probability densities. We compare the ML and 
Bayesian methods in an analysis of simulated data from Rhodamine 110 for various 
combinations of expected fluorescence emissions and background noise. 

Model Derivation 

We assume that the arrival times of the photons from a single fluorescent 
molecule are recorded in a finite interval [0, T] . Let 7 , be the arrival time of the i' h 
photon, where 0 < t, < T . We assume that the number of photons detected from the 
fluorescent species is a Poisson process with continuous intensity function X(t) = Ae ~ th , 

for t> 0 . The parameter t is the average fluorescent lifetime and A is the number of 
photons detected at time 0. Let A(t) = X(u)du . The quantity A(T) is the average 

fluorescence emission in the observation interval (0, T] . We assume that the number of 
background photons detected is a Poisson process with a continuous intensity function 
<7(0. t >0 . Forms of q(t) which have been reported include the constant function and a 

function composed of a linear combination of two Gaussian functions, and exponential 
function and a constant. Let Q(t)= ^q(u)du- We assume that q(t) is known from 

background photon measurements made on the detection system in the absence of the 
fluorescent molecules. We assume also that 



b -> oo 


Equation (1) is a technical assumption which ensures that the joint, conditional and 
marginal probability densities of the arrival times are well defined. Under the assumption 
that the fluorescence emission and background noise processes are independent, the 
number of photons arriving from both is a Poisson process with intensity function 
^(0 + < 7 ( 0 - 
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To derive the joint probability density of photon arrival times we consider the 
event that / photons arrive in (0, T] at times 0 < t x < t,<T. This probability density 

is defined by the events that no arrivals occur in the intervals (0,/|],(/. + a/ ( ,/ i+1 ] for 
/ = and (t,+At h T], and that exactly one arrival occurs in each of the intervals 

(/,,/, +A/ ; ] for / = By the definition of an inhomogeneous a Poisson process the 

following statements describe the probabilities of these events: 

Pr (No arrival in (0,/, ]) = exp{- K(u) + q(u)du} 

Pr(No arrival in (l, + A/, +1 ]) = exp{- f' + ' X(u) + q(u)du) (2) 

v + A/j 

, K/+A// tf,+A/ y 

Pr (One arrival in (t, , t, + At, ]) = I X(u) + q(n)du exp{- I X(u) + q(u)du} 


The intervals are nonoverlapping and thus, independent by the basic axioms of a Poisson 
process. The joint probability of the events is therefore, 


/('|t,-4)TTa/,= exp{- f'^iO + ^zO^rrexpl- f' + ' X(u) +q(u)du} exp{- f X(u) + q(u)du} 
/= 1 10 7 = , *1-1 +A, | 4/ +A(/ 

( 3 ) 

i — r rt ^ a/, r i ^ —L— 

X I 1 4 M«) + < 7 (H)rf«exp{-J \(u) + q(u)du}+ o(| jA/y) 

where t^(t h t 2 ,...,t,) T . Dividing both sides by J“[a^ and letting A/ ; -»0 for each i gives 

i=i 

the joint probability density of the arrival times 


/(*|t ,/0 = J^WVT + ^Tlexpf- P A.(w) + ^(«)^i/}exp{- f \(u)+q(u)du} 

i = I * 

I-\ . 

xn«P(-P ' H u ) + q(u)du} (4a) 
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= [*■(“)+ •?(«)<*<} 
/=1 ■° 


for 0 </,<>...,<//< T , where we used the fact that 


(4b) 


r 


+ A/, 


X(w) + q(u)du 


At; 




( 5 ) 


as A/j — > 0 , for all i. 

Equation (4) shows that the joint probability density of the photon arrivals may be 
represented without the need to bin the data. This derivation may be viewed as taking the 
number of bins equal to the number of detected photons, placing the i th photon arrival in a 
bin of width A t t , and letting the widths of all bins go to zero. From equation (4b) it 
appears that the /, ’s are a collection of independent observations on the interval (0,71, 
however, they are not. Their joint probability density is defined on r ' subject to the 
constraint that Ocq < T<ao. In particular, equation (4a) shows that the arrival 

times are a Markov process in that the conditional probability density of t, depends only 
on . That is, we may write 



where fOiUi-i) = (M*,-) + ?(*<)) exp{-f' X(u) + q(u)du) and t. is defined on the interval 

■v-i ' 

('/-l. 00 ) . The practical reason for observing that the arrival times obey a Markov process 
is that this representation of the joint arrival time density suggests a simple method of 



page 8: Brown, Zhang, McCollom; TR 94-03; Single Molecule Fluorescence Measurements 

simulating photon arrival times for an inhomogeneous Poisson signal plus noise model 
(Lewis and Schedler, 1 977). This algorithm is presented in the Appendix. 

Conditional Arrival Time Probability Density 

The observed photon arrival time data may also be analyzed in terms of the 
conditional arrival time probability density. That is, given I photon arrivals detected in 
(0, 7] , what are the probable locations of the arrivals in that interval? The number of 
photons arriving in (0,71, N{T) , is distributed as a Poisson random variable with parameter 
\{T) . From equation (4) and the properties of a Poisson process the conditional arrival 
time probability density is 

/(t|T, A, /) = Pr (0 < ( t < T and N(T) = 7} l Pr {N(T) = /} 

/ 

]“[[X(^) + ^(/,)]expH 0 r k(u) + q{u)du) 

= t (6) 

[ A(T) + QiT)]' exp{-[A(73 + Q(T)]} / II 

I 

_ t = i 

[A(T) + Q(T)]' 

In the case where the Poisson intensity parameter is constant and the observation interval 
is divided into k bins, equation (6) is a ^-nominal probability mass function. 


Statistical Analysis of Experimental Data 

Maximum Likelihood Approximation of the Fluorescent Lifetime and Energy Emission 
Probability Densities. 
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From the joint arrival time density we can estimate t and A by the method 
maximum likelihood. The log likelihood based on the joint arrival time density is 

i 

log / 01 t ,A) = Y logW/,) + q(tj))-[A(T) + Q(T)] (7) 

/=l 

and the ML estimates can be obtained numerically by finding the values of t and A 
which maximize Equation (7). Similarly, the log likelihood for the conditional arrival 
time density in equation (6) is 

i 

log f{t\ x,A,I) = £ logC^O, ) + q(t , )) - / log[ A(70 + Q(T)] (8) 

i=l 

+ log(/!). 


Equation (8) may also be used for ML estimation of t and A . 

Under the assumption that A is large we can derive approximate probability 
densities for x and A . Let = {x,k) T be the ML estimates of x and A respectively and 
define * = (x, A) . The observed information matrix G is the 2x2 matrix whose elements 

are c „ , * and It follows 

dr oidA dA z 

from the large sample theory of ML estimates that the approximate joint probability 
density of x and A is the Gaussian density defined as 


/iOM[0 = - 


■exp{-4(*-|i) 7 V '(at- jx)} 


2n\V\ 2 


( 9 ) 


where V = -G 1 and \V\ is the determinant of V . The marginal probability density of 
/(X 10 , is the Gaussian density with mean x and variance equal to \\ , , i.e. the /, / element 
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of V . To compute the approximate probability density of A(F) , we define the 
transformation 



( 10 ) 


and its Jacobian 


/ = 


A[ l-(l + ^)e- r/t ] T(I- e " m ) 


The determinant of J is ]/[= x(\-e~ Th ) . Applying the change of variables in equation 
(10) to /,(t, A\t) under the assumption that is A large, shows that the approximate joint 
probability density of t and A(T) is the Gaussian density defined as 


/ 2 (t,A(710 = — 

2n| S| 2 


( 11 ) 


where E =JVJ T . The approximate marginal probability density of A(f) is the Gaussian 
density with mean A(T) and variance I 2 2 ■ 


Bayesian Estimation of the Fluorescent Lifetime and Energy Emission Probability 
Densities. 

To conduct a Bayesian analysis of the lifetime and energy emission estimation 
problem we assume that knowledge about x and A known prior to the experiment can be 
summarized in terms of a prior probability density. Since the fluorescent lifetime is 
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independent of A , the number of photons observed at time 0, we write the joint prior 
density of x and A as the product of two locally uniform prior densities defined 
as/(x,/f) = /(t)/(/l) where, 


/(*) = 


0 


T, < T < x 2 

otherwise 


f(A) = 


(A 2 -Air' 
0 


A] < A < A2 
otherwise 


( 12 ) 


The values of and a 2 are defined from known properties of the particular single 

molecule experiment. Applying Bayes theorem equations (4) and (12) yields the joint 
posterior density of 




fit) 


(13) 


where f(t) = j J f(T)f(A)/(t\i, A)dAdx . 

T A 

Equation (12) shows that the locally uniform prior densities can be used to 
constrain the parameter valves to a region in the t - A plane which is physically plausible. 
Equation (13) shows that the uncertainty in the parameter defined by the Bayesian 
analysis is the product of the ranges defined by the prior densities and the uncertainty in 
or information about the parameters derived from the experiment, and summarized in the 
likelihood function f(t\\,A) . That is, the Bayesian analysis combines prior information 
with that collected in the experiment to determine the uncertainty in the model 
parameters. Because the prior density is uniform on the t -A plane and because the 
likelihood and log likelihood are equivalent summaries of the experimental information. 
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if the experiment contains a lot of information about the model parameters then the 
uncertainty in the model parameters defined by the ML and Bayesian analyses will agree. 

The posterior probability density of the fluorescent lifetime is obtained by 
integrating (13) with respect to A to yielding 

(14) 

A 

Since we are interested in A(T) instead of A we apply the change of variables defined by 
equation (10) and find that the joint posterior probability density of the fluorescent 
lifetime and fluorescent emission is 

hi t, A(r)|0 = -W,(t, " «" r/T )l') ( 1 5) 

I J I 

The marginal posterior density of A(7") is defined as 

/ 2 (A(r)|0= J/ 2 (t,A(r)IO* (16) 

Equations (14) and (16) define the uncertainty in the fluorescent lifetime and energy 
emission respectively in terms of probability densities without having to make large 
sample assumptions. The respective modes, medians and means of /i( T IO and / 2 (A(r)|0 
may be used as point estimates of x and A(T) . 


Simulation Study of Rhodamine 110 Single Molecule Fluorescence Detection 

We compared the ML and Bayesian methods in a study of simulated single 
molecule fluorescence experiments for Rhodamine 1 1 0 assuming an expected fluorescent 
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lifetime and detection interval of 4.2 and 10 nanoseconds respectively based on 
Tellinghuisen, et. al. We simulated three levels of expected background noise fraction: 
high = 0.70, moderate = 0.35 and none = 0 combined with four levels of expected 
fluorescence emission: low = 19.1, medium =191, moderate =952 and high =3812 
photons 10 nanoseconds. These expected fluorescence emission levels are the values of 
A(T) obtained with A set respectively to 5, 50, 250 and 1000 with t =4.2 and r=10. We 
assumed q(t) = k n , a constant function, whose value was determined for any given 
specification of the noise fraction and expected fluorescence by the relation 


MT) + x„t 

where nf is the expected background noise fraction and X n T = Q(T ) . 

For each combination of nf and A(7")we simulated 3 detection experiments (36 in 
total) using the algorithm defined in the Appendix. ML estimation was carried out by a 
combination quasi-Newton’s method and local search procedure and the Bayesian 
methods by rectangular integration procedures. For the Bayesian analysis we took 
x, = o and x 2 = 12 nanoseconds in the prior density for /(x) . Given previous reports on the 
fluorescent lifetime of Rhodamine 1 10, it is reasonable to assume that x would lie in this 
interval. For f{A) we chose four different representations for both A\ and A 2 depending on 
value of A used in the • simulations. These were 


A 

4 

a 2 

5 

0 

10 

50 

0 

100 

250 

150 

350 

1000 

750 

1250 



page 14: Brown, Zhang, McColIom; TR 94-03; Single Molecule Fluorescence Measurements 


Maximum Likelihood Analysis 

As nf decreased and ef increased, the precision in both the fluorescent lifetime and 
energy emission estimates increased (Table 1). The CV's of the lifetime and energy 
emission estimates decreased respectively from 58.03% and 46.8% ( nf = 0.35 and ef = 
19.1) to 2.68% and 1.61% («/= 0.0 and ef= 3812). In all the simulations the true values 
of t and A(T) are covered by the 0.95 confidence intervals based on the maximum 
likelihood estimates (mean ± 2sd). The widths of the 95% confidence intervals (4sd) 
fori ranged from 8.4 ns {nf = 0.35 and ef = 19.1) to 0.44 ns ( nf =0.0 and ef = 3812). 
The absolute widths of the 0.95 confidence intervals forA(7) increased with increasing 
values of A{T) however, the length of the interval as a percentage of A {T) decreased 
from 93.6% (nf = 0.35 and ef = 19.1) to 6.6% ( nf =0.0 and ef= 3812) , i.e. twice the 
CV's. Due to their low signal to noise ratio the simulations wi*th ef= 19.1 gave point 
estimates of t which agreed least with the intended value of 4.2 ns. The increase in 
estimation precision of xandA(r) with decrease in nf and increase in ef is shown in 
Figure 1. The estimated joint density for nf = 0.35 and e/= 19.1 has approximately 10% 
of its support on negative values of t and A(T) because the Gaussian approximation 
based on the ML estimates is not restricted to positive values (Figure 1A). Across all the 
simulations the approximate joint densities of t and A(T) have slight positive 
correlations (0.3 1 to 0.35) for nf = 0.70 (Figures 1 A, D and G) whereas the two estimates 
were uncorrelated for nf = 0.0 (Figures 1C, F and I). 


Bayesian Analysis 

For the Bayesian analysis as nf decreased and ef increased, the precision in both 
the fluorescent lifetime and energy emission estimates also increased (Tables 2 and 3 and 
Figure 2). The CV's of the lifetime and energy emission estimates decreased respectively 
from 55.03% and 41.51% (nf = 0.35 and ef= 19.1) to 2.69% and 1.61% (»/= 0.0 and ef= 
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3812). In all the simulations the true values of x and A(7)are covered by the 0.95 
credibility intervals, i.e. the Bayesian equivalent of the 0.95 confidence intervals. For 
ef = 952 and 3812 there was very good agreement between the joint probability densities 
ofx and A(7) and each of their marginal probability densities (Tables 1, 2 and 3). All the 
Bayesian probability densities for these values of ef appeared to be Gaussian as 
suggested by their plots and the fact that their skewness and kurtosis values are all close 
to zero. Similarly, for the simulations in which <?/= 191 and nf equal to 0.35 or 0.0, there 
was good agreement between the ML and Bayesian probability density estimates 
(Figures IE and 2F, and Figures 2E and 2F). 

For ef =19.1 for all values of nf and for ef =191 and nf = 0.70 the Bayesian 
estimates of the joint probability densities ofx and A {T) are non-Gaussian (Figures 2A, 
2B, 2D and 2G). The marginal probability densities for A(7) are nearly Gaussian for these 
four combinations of ef and nf and these densities are in very good agreement with their 
ML counterparts. The case in which the marginal density is least Gaussian is shown in 
Figure 3D. The non-Gaussian nature of the joint probability density stems mostly from 
the uncertainty in x as shown in Figures 2A-2D and indicated by the non-zero values for 
the skewness and kurtosis obtained in these simulations (Table 2). All four of the 
marginal probability densities of x have tails which are skewed to the right, i.e. positive 
skewness. In the simulations for ef =19.1 and nf = 0.70 and 0.35, the upper bound of 
these probability densities is defined by the prior density. The 0.95 credibility intervals 
for x in these two simulations extends from 0.14 to 11.4 ns and from 2.3 to 11.7 ns. 
These findings suggest that under these experimental conditions the data contain only 
minimal information about the molecule's fluorescent lifetime. 
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Discussion 

Our approach provides an exact probability model of the photon arrival times 
observed with background noise and makes explicit the relation between the counting 
process, arrival time and conditional arrival time probability densities. This model avoids 
the need to bin data and extends the analysis provided [] and [] by making possible ML 
and Bayesian analyses which consider the background noise in the estimation of the 
excited species lifetime and energy emission. Both procedures measure the uncertainty in 
these quantities based on data from the photons collected from single molecule in a single 
experiment. 

Our ML analysis gives explicit formulae for the change in precision of the 
lifetime estimates as a function of the signal intensity, background noise and detection 
interval. ML analysis of experimental data provides the most efficient use of 
experimental information and as such, is the preferred approach provided the ML 
computations are numerically tractable the experiment is highly informative about the 
model parameters relative to the known prior to the experiment. Inferences about 
uncertainty in model parameters can be made from ML theory as long as the appropriate 
asymptomatic assumptions are satisfied. Contrary to the suggestion by Tellinghuisen et 
al., our results show that ML estimation for single molecule lifetime and energy emission 
and is highly tractable procedure which can be carried out efficiently using a quasi- 
Newton’s methods. Hall and Selinger showed that when the observation interval is finite 
the large sample properties of the ML estimate of t and A are not of order / but of order 
A . Our ML analysis extends their result by showing that the result also hold when the 
data consist of a Poisson signal plus Poisson background noise model observed in finite 
time. 

While the ML and Bayesian procedures will give similar parameter and 
uncertainty estimates when the expected signal intensity is large, our simulation study 
shows that when the number of detected photons is small, the Bayesian procedures 
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provide a more reliable uncertainty assessment of uncertainty. The Bayesian procedure 
makes it possible to combine prior information about the most probable range to x and A 
prior to the experiment with the information on these parameters summarized in the 
likelihood. In particular, the Bayesian procedure trades the large sample assumptions 
required by the ML error analysis for the more realistic one that the fluorescent lifetime 
and energy emission— where the latter is implicitly defined by x and A —lie within given 
intervals. As we have shown, this interval may be specified based on the known 
properties of the experimental paradigm. The Bayesian procedure suggests therefore, a 
preferred alternative when the expected signal intensity is low and as a consequence, the 
large sample assumptions needed to justify the ML analysis cannot be satisfied. 

While either the Bayesian or ML procedures works well for analyzing 
fluorescence emissions, the Bayesian methods provide more realistic measures of 
uncertainty in the fluorescent lifetimes for any combination of background noise and 
fluorescence emission. The Bayesian methods should provide more realistic assessments 
of uncertainty in fluorescent lifetime estimates in current single molecule flow stream 
experiments where the expected fluorescence emission can be in the low to medium 
range. Both the ML and Bayesian algorithms can be automated for real-time 
applications. 
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Appendix 

To derive the algorithm we note that the conditional cumulative distribution function 
(CDF) of t, given is 


^0/K-l) = 1 -exp X (u) + q(u)du} 


(A.l) 


Equation (1) ensures that FO/hy-i) is a proper CDF defined on ») . Sequential 
application of the inverse CDF algorithm for Monte Carlo sampling gives the following 
algorithm: 

STEP 0: Set t 0 = 0 . 

STEP 1 : Given draw r t a uniform random number on (0,1) . 

STEP 2: Compute F~\r,) = r from (A.l). 

STEP 3: If i' < j then set t, = , = t, and go to STEP I, otherwise stop. 

The inversion in STEP 2 is computed numerically. 
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TABLE 1 . 


MAXIMUM LIKELIHOOD ESTIMATION SUMMARY: 


FLUORESCENT LIFETIME 
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ENERGY EMISSION 


CVxl00% 

MEAN 

SD 

CVxl00% 

58.03 

16.60 

7.77 

46.81 

44.97 

19.60 

5.49 

27.99 

27.44 

22.10 

4.71 

21.32 

20.82 

207.25 

25.02 

12.07 

15.29 

221.25 

17.98 

8.12 

11.22 

182.75 

13.55 

7.41 

13.25 

876.00 

55.79 

6.37 

6.88 

964.00 

38.19 

3.96 

5.31 

908.00 

30.13 

3.30 

2.87 

3950.00 

63.50 

1.61 

2.93 

3880.00 

65.78 

1.62 

2.68 

3870.00 

62.19 

1.61 



TABLE 2 . 

BAYESIAN ANALYSIS SUMMARY: FLUORESCENT LIFETIMES 
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tThe skewness of a probability density of a random variable X is the normalized third moment defined as £(Y-n) 3 /a 3 where n is the mean and 
is the variance. 

ttThe kurtosis of a probability density of a random variable X is the normalized fourth moment defined as E(X - n ) 4 / a 4 - 3 . 
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